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Abstract. Collisions of cosmic rays in the atmosphere may produce heavy hadrons of very 
high energy. The decay length of a B meson of energy above 10 7 GeV is larger than 1 
km, implying that such a particle tends to interact in the air before it decays. We show 
that the fraction of energy deposited in these interactions is much smaller than in proton 
and pion collisions. We parameterize their elasticity and determine the average number of 
interactions and the atmospheric depth at the decay point for different initial energies. We 
find that the profile of a 3 x 10 9 GeV bottom shower may be very different from the profile of 
a proton shower of the same energy, defining either a very deep maximum, or two maxima, 
or other features that cannot be parameterized with a single Gaisser-Hillas function. Finally, 
we discuss under what conditions a bottom hadron inside the parent air shower may provide 
observable effects. 
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1 Introduction 

Cosmic rays produce collisions in the upper atmosphere of energy well above the ones observed 
in particle colliders. Our knowledge of the hadron properties must then be extrapolated to 
such energies, a procedure that introduces some important uncertainties. These uncertainties 
are not only related to the appearance of new particles and interactions, but also to the 
transition to a regime where the properties of the standard particles could be substantially 
different. In particular, heavy hadrons that decay through weak interactions have typical 
decay lengths Aq 6C = 0.1-0.5 mm. As their energy grows, Lorentz dilation will eventually 
make Aq gc longer than the interaction length in the air, implying that they tend to collide 
before decaying. Obviously, we do not know how a charmed or a bottom hadron behaves 
inside a calorimeter, since in colliders they decay well before they can reach it. Although 
heavy-hadron interactions are of no interest there, they are processes that might occur in 
extensive air showers and could introduce unexpected effects. The production of standard 
heavy-hadrons in the atmosphere has been extensively studied in the literature (see [1] for a 
review), as its prompt decay is an important source of muons of very high energy (see also 
[2, 3]), however, their propagation has been basically ignored. 

In a high-energy collision with a target at rest, a proton or a pion will break into several 
pieces, with each piece taking a similar amount of energy. The elasticity z (fraction of energy 
taken by the leading hadron) and the multiplicity in such collisions will determine how the 
hadronic shower develops through the atmosphere, which is equivalent to 10 meters of water 
if crossed vertically. It is clear, however, that if the collision involves a charm or, especially, 
a bottom hadron the situation will be very different, as the piece carrying the heavy quark 
will take most of the initial energy after the collision. Therefore, qualitatively one expects 
lower inelasticities and, as a consequence, a different profile of deposited energy. 

In this paper we will explore the propagation of b hadrons in the atmosphere. In Section 

2 we will follow the method applied for charmed hadrons in [4] to parameterize the elasticity 
of their collisions with air nuclei. In Section 3 we use CORSIKA [5] to find the profile of 
bottom showers of different energies. In Section 4 we discuss under what circumstances a 
bottom component inside an extensive air shower may give an observable effect. Section 5 is 
devoted to the summary and discussion. 
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2 Inelasticity in heavy-hadron collisions 

In this section we will modify PYTHIA [6] to simulate the collisions of a hadron H (contain- 
ing a b quark, H = B) with matter. PYTHIA describes hadronic collisions of pions and 
protons, distinguishing two types of interactions: diffractive processes, where the two hadrons 
(as a whole) exchange momentum through pomerons, and non- diffractive or partonic ones, 
where gluons are exchanged between the partons in the colliding hadrons. The latter type 
includes both soft collisions of q 2 < 1 GeV and the hard ones in deep inelastic scatterings. 

In PYTHIA non-diffractive processes dominate the inelastic cross section. For example, 
the cross section of a 10 9 GeV proton or pion with a proton at rest is 

o- pp = 94.1mb; = 64.8 mb , (2.1) 

where diffractive processes contribute just 30% in pp collisions and 26% in 7rp interactions. 
Therefore, we will study diffractive and partonic processes separately, then we will combine 
them to define a generic Hp scattering, and finally we will discuss the collision with an air 
nucleus. 

2.1 Diffractive Hp collisions 

Let us first describe the approach to simulate a diffractive Hp collision with PYTHIA from 
the analogous light-hadron process. A hadron H with mass mjj contains a heavy core of mass 
mi, ~ 4.7 GeV. One may think of a proton as three clouds of mass m = 0.3 GeV associated 
to the three constituent quarks. Within that picture, a A& baryon will consist also of three 
similar clouds, but with one of them having an additional electroweak core and a total mass 
mb + m = 5.0 GeV. In a diffractive process, however, this heavy core will be invisible (the 
proton and the A& clouds will look identical), since the interaction has a q 2 much smaller 
than 1 GeV 2 and cannot resolve it. Therefore, the momentum exchanged through pomerons 
or other non-perturbative dynamics with the target should not depend on the electroweak 
core in H. 

More precisely, the light degrees of freedom in H carry just a fraction w = [rein — 
fnb)l m H of the hadron energy E. In a diffractive scattering, H will be seen by the target 
nucleon like a light hadron of energy wE. Therefore, to estimate the momentum absorbed 
by H we will just simulate with PYTHIA the collision of a proton (for H = Af,) or a pion (for 
H = B) of energy wE with the nucleon and assume that q^ is the same when the incident 
particle is the heavy hadron. 

This momentum is all we need to simulate how H evolves after the diffractive collision 
(see [4] for details). Once H absorbs q^, it becomes a system of mass M* , the critical 
parameter in the collision. If M* < mu + 1 GeV the process is quasi-elastic and Hm? will 
just decay into two bodies (e.g., H + rj). For larger values of M* the system is treated by 
PYTHIA [7] like a string with the quantum numbers of H. When H is a baryon the string 
may be stretched between a quark and a diquark or between a quark, a gluon and a diquark, 
whereas for a diffractive meson the string connects a quark and an antiquark or a quark, a 
gluon and an antiquark. 

In Fig. 1 we plot the distribution of the fraction z of energy taken by the bottom hadron 
after a diffractive collision of a 10 9 GeV A^ or a B with a proton at rest. The average values 
are (z) = 0.88 and (z) = 0.90, respectively, which imply an inelasticity K = 1 — (z): 

Kfl =0.12; rf = 0.10. (2.2) 
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Figure 1. Elasticity (z) in Af,p (solid) and Bp (dashed) diffractivc collisions. 



In comparison, we note that for proton and charged-pion collisions, the inelasticity obtained 
also with PYTHIA is substantially higher: 

K$ = 0.40 ; Kfl = 0.23 . (2.3) 

2.2 Partonic Hp collisions 

For non-diffractive collisions, we model H as a system with the same parton content as the 
corresponding proton or pion, but substituting a valence up quark (uq) for a bottom quark. 
As in diffractive processes, we will associate a hadron H of energy £ to a light hadron of 
energy (mjj — m^)/m,H E. If uq carries a fraction x of the proton or the pion momentum, we 
will change it for a b with momentum fraction xi, 

x b = ^ + ^^x. (2.4) 

In this way the excess of energy in H is carried entirely by the bottom quark, whereas the 
light partons in both hadrons (H and p or ir) carry exactly the same amount of energy. 

We will then distinguish two types of Hp partonic collisions: those where the bottom 
is an spectator (i.e., it is a light parton in H which hits a parton in the target proton), and 
processes where the b quark itself interacts. For the first case we first just simulate with 
PYTHIA the parton process using a light hadron, later we then substitute the spectator uq 
for the bottom quark. Intrinsic bottom interactions, on the other hand, have a much smaller 
cross section, 

^L = 0-8mb, (2.5) 
than the processes with an spectator b quark, 

^pec = 56.2mb; a^ pec = 38.7 mb , (2.6) 
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Figure 2. Elasticity (z) in A;,p (solid) and Bp (dashed) non-diffractive collisions. 

but they imply collisions of higher inelasticity. 

In Fig. 2 we plot the distribution of energy taken by the bottom hadron after a A(,p or 
a Bp non-diffractive collision. The average inelasticity is 

K n-dif = Q 22 . K n-d lf = Q 21 ^ (2 7) 

For comparison for protons and pions of the same energy is: 

= 0.66 ; K%p dif = 0.77 . (2.8) 
2.3 Inelasticity in Hp collisions 

We need to combine both types of collisions. We find that the total (inelastic) Hp cross 
section for a projectile of energy 10 9 GeV is 

^ P = 4} P + * 4 / = 82.1mb; 

^ P = ^/ P + ^/ = 54.4mb, (2.9) 

where diffr active processes contribute 30% in A;,p collisions and 28% in Bp interactions. The 
cross sections for proton and pion collisions of the same energy are 14% and 19% larger, 
respectively. 

Using the relative frequency of these processes to find the average inelasticity, we obtain 
(see Fig. 3) 

K AbP = 0.19; K Bp = 0.18, (2.10) 
which is substantially lower than the inelasticity in proton and pion collisions, 

K pp = 0.59 ; K np = 0.70 . (2.11) 
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Figure 3. Elasticity (z) in (Ab,B)-p and (p, 7r)-p collisions. 
2.4 i/ air collisions 

The total cross section for the collision of H with an atomic nucleus of mass number A can 
be approximated as 

a HA ^A 2 / 3 a H ?, (2.12) 

where the factor of A 2 / 3 takes into account the screening between the nucleons inside the 
nucleus. For an averaged atmospheric nucleus of A = 14.6 this implies 

a A b air = 49Q mb . a Bair = 325 mb _ ^.13) 
The associated interaction length X? t = m a i r / a H air in the atmosphere is therefore 

At 6 i = 49gcm- 2 ; \f nt = 74 g cm" 2 , (2.14) 

which is a 14% and a 19% longer than those of a pion and a proton of the same energy, 
respectively. 

To deduce the energy and the species of the b hadron after the collision we follow [4] 
and distinguish between peripheral and central collisions. We assume that the spectrum in 
peripheral processes coincides with the one in Hp collisions, whereas the average inelasticity 
in central processes is the typical in a partonic scattering increased by 10%. In addition, we 
will take equal frequency for both types of processes. The average inelasticity can then be 
estimated as 

K Hatr » /rac!2 Kg*, + \ & « Q K Hp + \ 1.1 . (2.15) 

At E = 10 9 GeV, we obtain 

K\ b air ~ 0-21 ; K Bair ^0.20. (2.16) 
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Figure 4. Elasticity (z) in Af,-air (thick solid) and B-air (thick dashed) collisions. We include in 
both cases the fit with the exponential in Eq. 2.18. 

For proton and pion collisions the same prescription gives an inelasticity 

K pair « 0.66 ; K nair « 0.78 . (2.17) 

This 12% increase in K when going from a proton to a nucleus target compares well with 
the results obtained by other authors [8]. 

In Fig. 4 we plot (thick lines) the final distribution z of the fraction of energy taken by 
the heavy hadron after the collision of a A& or a B with an air nucleus. We also plot (thin 
lines) our parameterization of the distribution as an exponential between and 1: 

f(z) = l e -»(—°) 2 , (2.18) 

where N normalizes the distribution to 1. We obtain a = 31.3, zq = 0.799 and N = 0.30 for 
an incident A& baryon and a = 37.0, zq = 0.800 and N = 0.28 for an incident B meson. 

Finally, we plot in Fig. 5 the frequency of the different hadron species after the collision 
of a 10 9 GeV B meson with a proton at rest for different values of the elasticity z. We obtain 
a B° 42% of the times, B~~ mesons appear after 41% of the collisions, whereas the frequency 
of B® mesons and A^ baryons is 9% and 7%, respectively. In A& collisions the approximate 
frequency of (A&, B°, B~, B°) provided by PYTHIA is (68%, 14%, 14%, 4%). 

3 Propagation in the atmosphere 

The propagation of a long-lived b hadron differs from that of a proton or a pion, due to 
its larger elasticity in collisions with air nuclei and its longer interaction length. Using the 
Monte Carlo code CORSIKA we analyze the energy deposition of a b hadron produced in 
the upper atmosphere, detached from the development of its parent shower. We focus on B° 
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Figure 5. Frequency of the different 6-hadron species after the collision of a 10 9 GeV B° meson with 
a proton at rest. 

mesons, as baryons tend to convert into mesons after a few hadronic interactions. The B° 
meson may interact 1 (according to the cross section and the inelasticity described in Section 
2) or decay 2 as it propagates, producing secondary showers at different atmospheric depths. 

As an example, we will take a B° meson produced at a depth of Xq = 100 g cm -2 
with a zenith inclination 9 = 60°. We consider different values of the initial energy of the 
B° meson, Eb- For values of Eb of the order of 10 8 GeV, the B° mesons decay after few 
interactions and therefore look like a proton shower. When larger energies O(10 9 GeV) are 
considered, the differences between B° and proton initiated showers become more evident. 
To stand out the most relevant features of the showers produced by B° particles, we consider 
in this section average proton shower profiles. The next section, devoted to the detectability 
of bottom hadrons produced inside extensive air showers, will deal with individual fully- 
simulated showers to take into account the effect that shower fluctuations introduce in our 
analysis. 

In Fig. 6 we plot the averaged profile of 100 B° initiated showers (solid) and of the same 
number of proton showers (dashed). The initial energy for B° and protons is 3x 10 9 GeV. The 
average B° decay takes place at (Xf) = 1190 ± 340 g cm -2 , after (Nj) = 16 ± 4 interactions 
with air nuclei. These B° showers are different from a proton shower of the same energy 
in two important aspects. First, the energy deposition rate in the atmosphere is slower, 
and the shower maximum is reached later, around 1000 g cm -2 . Among the 100 showers 
generated, 52 reach their maximum beyond 1000 g cm -2 . Also, secondary showers of energy 
10 7 -10 8 GeV start deep in the atmosphere (as a result of the decay or the collisions of the 
heavy hadron), increasing the number of electromagnetic particles that reach the ground. 

If it deposits an energy AE = (1 — z)E in a collision, we will simulate there a CORSIKA pion shower of 
the same energy AE. 

2 If it decays at a depth X we treat the decay using PYTHIA and inject the products at that depth using 
CORSIKA. 




-7- 



xlO 6 



100 






n — j x ±u ucv 

\ 


80 




P/' 

• 




d 




• 
* 

* 
* 


* 


ni 60 




* 
* 




>< 




■ 
• 












uj 40 












• / 


» \. 

* X 
» \^ 


20 




• 1 


» \» 
* 

V 







* M 


, . . . , . . . -r-.. 



400 800 1200 1600 
depth [g cm -2 ] 

Figure 6. Average profile of a p (dashed) and a B (solid) shower for an initial energy of 3 x 10 9 GeV. 




Figure 7. Anomalous profile of a 3 x 10 9 GeV B° shower (histogram). The solid line corresponds to 
a fit using a combination of two Gaisscr-Hillas functions. 

The profile of these showers, in general, cannot be well described by a single Gaisser-Hillas 
function. 

In Fig. 7 we plot as an illustration an anomalous profile found among these 100 showers, 
namely a longitudinal development featuring a double peak topology. We obtain this feature 
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Figure 8. Profile of a 2 x 10 10 GeV proton shower that has produced a 3 x 10 9 GeV B meson at 
X Q = 100 g cm" 2 . 

in 15 out of 100 simulations. In the same figure we plot (solid) a fit to this shower using the 
sum of two Gaisser-Hillas. 

4 Detectability of B mesons inside extensive air showers 

B mesons are not primary particles that may start an air shower. They will be produced 
in the collision of a primary cosmic ray or an energetic leading hadron in the shower with 
an air nucleus. Therefore, the B° shower described in the previous section will always be a 
component inside a parent shower. Its observability will then depend critically on the fraction 
Xb of the energy that the B meson carries when it is produced. Values x^ > 0.1 are interesting. 
In Fig. 8, we plot for the sake of the argument an extreme example where a CORSIKA- 
simulated proton primary with energy E p = 2 x 10 10 GeV and zenith angle 9 = 60° produces 
a very energetic B° meson in the upper atmosphere, carrying Eb = 3 x 10 9 GeV (xf, = 0.15). 
The rest of the secondary particles produce a shower with energy E p > = 1.7 x 10 10 GeV that 
reaches a maximum around 800 g cm" 2 . At ground level the electromagnetic component of 
this shower has been almost completely absorbed. The B° meson, produced at Xq = 100 g 
cm" 2 and decaying at X = 1748 g cm" 2 (depth along the proton incoming direction), has a 
deep maximum and a rich electromagnetic component at the ground. We observe a long tail 
in the shower profile, with dE/dx values close to the ground that are anomalous in a regular 
proton shower. 

Based on the model described in this work, the simulations show that the inclusion of 
b hadrons in the shower does not produce a striking deformation of the shower profile, but 
rather some features rarely found in the development of proton showers. Even so, these show- 
ers will be very often similar to those of protons. Thus, the search for a bottom component 
should be done on the basis of large statistics, trying to determine which are the observ- 
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ables that provide the maximal separation between proton showers and proton showers with 
bottom production. 

To assess the observability of b quarks inside inclined extensive air showers, we have 
produced several samples of fully simulated extensive air showers. These CORSIKA showers 
naturally include fluctuations in the longitudinal development of the cascades, which can 
mimic the effects produced by a bottom component. All in all, we have fully simulated 
several thousands of events for background and signal componentes: 2000 proton showers for 
each of the following energies logio(E/GeV)={10, 10.25, 10.5, 10.75} and fixed zenith angle 
of = 60° stand for our initial background sample. Per each of the former energies, another 
2000 proton showers with a bottom component constitute the signal. 

In our simulation, heavy quarks are produced only in the first interaction of the incoming 
proton with an air nucleus, and not in subsequent interactions of the leading or secondary 
hadrons with the air: their effect in the shower profile will be small due to their smaller 
energy. We have taken the cross section for bottom production in p-air collisions in [9] (they 
use a dipole picture in the Color Glass Condensate formalism, including saturation effects). 
It is found, for example, that at 10 10 GeV this cross section is around 25 mb, with the bottom 
component carrying more than 1% of the initial energy (i.e. x b i > 0.01) in a 25% of the 
collisions. Therefore, in these showers there would be around 1% probability to obtain a 
bottom component of energy above 10 s GeV in the first proton interaction. In our analysis 
we will take a total set of n showers, where Sq = 0.01 n showers have produced bottom 
hadrons (with the energy distribution in [9]) and -Bo = 0.99 n are regular proton showers. 
The propagation of the bottom component in these showers is performed according to the 
results obtained in Section 3. 

A thorough examination on an event by event basis of the simulated samples revealed 
a set of observables, related to the longitudinal development of the electromagnetic part of 
the shower, that help discriminating signal from background. The most relevant are: 

• The logarithm of the number of electromagnetic particles reaching the ground. 

• The difference between the depths beyond X max at which the shower is at 40% and 
35% of the shower maximum. 

• The difference between the depths before X max at which the shower is at 70% and 75% 
of the shower maximum. 

• The full width at 5% of the profile maximum. 

• The central moments of 3 rd and 5 th order of the shower profile. 

• The depth beyond X max at which the curvature of the Gaisser-Hillas fit to the profile 
changes. 

To enhance the separation power between signal and background, we trained a Fisher dis- 
criminant using the former observables [10]. Applying this discriminant and maximizing 
the value of S/V S + B, where S and B are the number of selected signal and background 
events, we obtain a signal selection efficiency E^ el = 0.165 and a background rejection effi- 
ciency E r ^ = 0.992. Thus, the number of selected signal and background events would be 
S = 0.00165 n and B = 0.00792 n, respectively. This means that, according to our model, to 
get a 3<T evidence for production of b hadrons in extensive air showers (with energy above 
~ 10 10 GeV and zenith angle between 55° and 65°) would require to collect around 3 x 10 4 
extensive air showers. 
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5 Summary and discussion 

Hadrons containing a b quark become long-lived at energies above 10 GeV, and they tend 
to collide in the air before decaying. We have shown that the inelasticity in those collisions is 
almost four times smaller than in light-hadron scatterings. As a consequence, these hadrons 
can cross a large fraction of atmosphere keeping a significant fraction of energy. For example, 
a 3 x 10 9 GeV B meson has an average of 16 interactions before decaying at Xf ~ 1200 g 
cm" 2 . The final energy of the meson is around 10 7 ' 5 GeV, and it will be all deposited close to 
the ground. We find that the profile of these showers frequently presents a very deep shower 
maximum and occasionally a double maximum. 

Although b hadrons are rarely produced in extensive air showers, we have also discussed 
whether they can provide a significant deviation once included in the shower development. 
We have shown that if the heavy meson decay occurs deep enough in the atmosphere the 
profile could reveal such a feature. A multivariate analysis of observables defined on the 
shower electromagnetic profile shows that, according to our model, a very large number of 
events need to be collected in order to get an evidence of b hadron production in ultra-high 
energy cosmic rays. 

Despite this result, we think that the production and propagation of heavy hadrons 
should be included in Monte Carlo codes like AIRES [11] or CORSIKA [5], which simulate 
extensive air showers of energy up to 10 11 GeV. Such an effort seems necessary to establish 
on a solid ground the expected frequency of events with late energy deposition and the best 
strategy in the search for observable effects. 
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